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ABSTRACT 



We present a self-consistent nonparametric model of the local cosmic ve- 
locity field derived from the distribution of IRAS galaxies in the PSCz red- 
O ' shift survey. The survey was analyzed using two independent methods, both 

' based on the assumptions of gravitational instability and linear biasing. The 

■ two methods, which give very similar results, were tested and calibrated on 

mock PSCz catalogues constructed from cosmological N-body simulations. 
The denser sampling provided by the PSCz survey compared to previous IRAS 
5— i ' galaxy surveys allows an improved reconstruction of the density and velocity 

fields out to large distances. The most striking feature of the model velocity 
field is a coherent large-scale streaming motion along the baseline connecting 
Perseus-Pisces, the Local Supercluster, the Great Attractor, and the Shapley 
Concentration. We find no evidence for back-infall onto the Great Attrac- 
tor. Instead, material behind and around the Great Attractor is inferred to 
be streaming towards the Shapley Concentration, aided by the compressional 
push of two large nearby underdensities. The PSCz model velocities compare 
well with those predicted from the 1.2Jy redshift survey of IRAS galaxies and, 
perhaps surprisingly, with those predicted from the distribution of AbcU/ACO 
clusters, out to 140 h^^ Mpc. Comparison of the real-space density fields (or, 
alternatively, the peculiar velocity fields) inferred from the PSCz and cluster 
catalogues gives a relative (linear) bias parameter between clusters and IRAS 
galaxies of be = 4.4 ± 0.6. Finally, we implement a likelihood analysis that 
uses all the available information on peculiar velocities in our local universe to 
estimate f3 = Qq-^ /b = O.G^Q^g (l"""); where b is the bias parameter for IRAS 
galaxies. 
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1 INTRODUCTION 

The art of modelling the cosmic velocity field, which originates from the desire to interpret observed deviations from 
a uniform Hubble expansion, has developed rapidly over the past few years. There are two main reasons for this. One 
is an increase in the quantity and quality of measured galaxy peculiar velocities. The other is the advent of nearly 
all-sky, flux-limited, redshift surveys that allow self-consistent theoretical predictions to be made with the requisite 
accuracy. 

Although other possibilities have been proposed (e.g. Babul et al. 1994), the gravitational instability theory 
(Peebles 1980) has proven to be the most successful theoretical framework in which to interpret peculiar velocity 
data in relation to inhomogeneities in the mass distribution. Early attempts to account for observed velocities within 
this general framework were rather simplistic due to incomplete knowledge of the distribution of galaxies in the local 
universe. Thus, simple parametric models were developed to describe cosmic flows in terms of infall onto one or more 
spherical overdensities such as the Virgo Cluster, the Great Attractor or the Perseus Pisces supercluster (Davis and 
Peebles 1983, Lynden-Bell et al. 1987, Peebles 1988, Han and Mould 1990). The situation changed dramatically when 
statistically complete, nearly all-sky, redshift catalogues of galaxies were constructed since it then became possible 
to predict peculiar velocities directly, assuming that luminous objects trace the underlying density fleld in some 
fashion. Since then, several modelling procedures have been developed which are generally based on the simplifying 
assumption that the gravitating mass is distributed just like the tracer objects (galaxies or clusters), although the 
relative amplitude of the deviations from uniformity, usually called the bias, is taken to be a free parameter. In 
addition to this "linear bias model," current methods also assume gravitational instability in the linear or mildly 
nonlinear regime (e.g. Yahil et al. 1991, Kaiser et al. 1991, Nusser and Davis 1994, Fisher et al. 1995b, Sigad et al. 
1998), and so require smoothing over scales where non-linear effects are important. This necessitates the additional 
assumption that smoothing a distribution that has evolved to a nonlinear state gives a result that can be modeled 
by linear or quasilinear evolution from smoothed initial conditions. 

Because of their large sky coverage, the most extensively used redshift surveys are those based on the "Point 
Source Catalogue" (PSC) of IRAS galaxies (e.g. the 1.9 Jy survey of Strauss et al. 1990, the upgraded 1.2 Jy catalogue 
of Fisher et al. 1995a, or the deeper but sparser QDOT survey of Rowan-Robinson et al. 1990). Other catalogues 
containing different kinds of objects such as optically selected galaxies (e.g. Shaya, TuUy and Pierce 1990, Hudson 
1994, Baker et al. 1998) or clusters of galaxies (Scaramella 1995 Branchini and Plionis 1996) have also been used to 
produce model velocity flelds on scales up to 200 h-^ MpcR. 

Comparison of a model velocity field with directly measured peculiar velocities provides a means to constrain the 
density parameter, f2o, while also testing the gravitational instability hypothesis and the assumed biasing scheme. 
A successful model for the peculiar velocity field may be used to recover the distribution of galaxies or clusters in 
real space, directly from their measured redshifts. This, in turn, allows investigation of the statistical properties of 
the objects' distribution, free from the effects of redshift space distortions. Overall, model velocity fields constructed 
using different surveys have proved to be remarkably consistent with one another (Freudling, da Costa and Pellegrini 
1994, Baker et al. 1998) and have succeeded in reproducing most of the characteristics exhibited by maps made 
directly from observed peculiar velocities. There are, however, two notable exceptions: the large bulk motion on very 
large scales claimed by Lauer and Postman (1994) and a dipolar coherence in the residuals between the velocity 
field obtained from the Mark III catalogue (Willick et al. 1997) and the IRAS 1.2Jy gravity field on scales of ~50 
ft-^ Mpc (Davis, Nusser and Willick 1996). 

Here we present a new nonparametric model of the cosmic velocity field based on the recently completed PSCz 
survey of IRAS galaxies. This is the last of the nearly all-sky redshift surveys based on the IRAS catalogue and 
supersedes both the 1.2Jy and the QDOT catalogues which it contains. The denser sampling and lower fiux limit 
of the PSCz survey allows us to model the peculiar velocity field on larger scales than was possible before without 
excessive contamination from shot noise. 

* Throughout this paper we write Hubble's constant as = lOO/i km s~^ Mpc~^ 
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The outline of this paper is as follows. In §2 we describe the PSCz dataset as well as two other redshift catalogues 
that we use to construct independent model velocity fields. Two methods for generating the PSC2: peculiar velocity 
fields have been implemented in order to keep track of systematic errors. These are presented in §3, together with 
a detailed error analysis. A cosmographic tour is performed in §4, along with a consistency check between the two 
PSCz velocity models and a comparison between the PSC2; and 1.2Jy model velocity fields. In §5 we take advantage 
of the large depth of the PSCz survey to compare the gravity field derived from it with the one derived from the 
distribution of Abell/ACO clusters. An estimate of the parameter /3 ~ 0,o'^/b, where b is the linear bias parameter 
of IRAS galaxies, is obtained in §6 by comparing observed and predicted bulk velocity vectors. In §7 we discuss our 
results further and summarise our main conclusions. 



2 THE DATASETS 

The main datasot used in this work is the recently completed PSCz redshift survey described in detail by Saunders 
et al. (1998a). The catalogue contains some 15,500 IRAS PSC galaxies with 60 /xm flux, /eo > 0.6. The construction 
of the catalogue is described in detail in Saunders et al. (1998a) and so we only summarize here its main features. 
The selection of potential galaxies from the PSC emphasised completeness at the expense of contamination by cirrus, 
stars, AGN, etc. Because the completeness of the PSC is not guaranteed to 0.6Jy in 2HC0N areas (Beichman et al. 
1988, PSC explanatory supplement), wc added to the PSC any IHCON detections from the Point Source Catalogue 
Reject file with FSS counterparts (Moshir ctal 1989). We then selected sources with /loo < 4feo (to eliminate the 
vast majority of cirrus sources) and /so > 0.6/25 (to eliminate the vast majority of stars). Residual contamination 
was eliminated by a combination of optical identification on Schmidt plates, examination of the raw IRAS addscan 
profile, information from the Simbad database and, where still unclear, optical spectroscopy. Large galaxies can have 
their flux underestimated by the IRAS beam, so wc constructed an all-sky list of optically-selected galaxies with 
extinction-corrected diameters larger than 2.25'. The PSC detections for these galaxies were flagged for deletion. The 
optical catalogue was addscanned by IPAC, and we used software provided by Amos Yahil to derive fluxes from the 
addscan profiles, and where the flux was sufficient, added them to the catalogue. 

For our purposes, one of the most important properties of the PSCz catalogue is its large sky coverage. The only 
excluded regions are two thin strips in ecliptic longitude that were not observed by IRAS, the Magellanic clouds, and 
the area in the galactic plane where the B-band extinction, Ab, exceeds 2 magnitudes. The extinction is estimated 
from the 100 //m background, and includes a correction for the estimated temperature gradient in the Milky Way 
(Saunders et al. 1998a). Overall, the PSC2 catalogue covers ~ 84% of the sky. Although we will occasionally consider 
galaxies with recession velocity as large as 30000 kms~^, for the most part we will restrict our analysis to the PSCz 
subsample of 11206 galaxies within 20000 kms~'^in the Local Group (LG) frame. The distribution on the sky of PSCz 
galaxies in shown in galactic coordinates (AitofF projection) in the upper part of Fig. 1. 

For comparison purposes, we have also applied our analysis to the similar, but shallower, IRAS 1.2Jy redshift 
survey (Fisher et al. 1995a). Galaxies in this catalogue were also selected from the IRAS PSC but using a higher 
(ADDscan) flux limit, /eo > 1.2 Jy, and different criteria for minimizing contamination by galactic cirrus. This 
catalogue has a somewhat larger sky coverage of ^ 88%. Here we will use the 4939 IRAS 1.2 Jy gala^xies within 20000 
kms~^of the Local Group. 

Finally, in an attempt to extend our analysis, we use a completely different set of mass tracers consisting of a 
volume-limited sample of optically-selected galaxy clusters extracted from the Abell (1958) and the Abell, Corwin 
& Olowin (1989; hereafter ACQ) catalogues. These were combined into a single homogeneous catalogue using the 
clusters in common between the two samples, according to the prescription described by Branchini and Plionis (1996; 
hereafter BP96). The sample has a limiting depth of 250 Mpcand contains 493 clusters of richness class R > 
at |6| > 13°, and mio < 17, where mio is the magnitude of the tenth brightest galaxy in the cluster. 

The number of galaxies in a flux limited sample decreases with distance, as may be seen in the lower plot of 
Fig. 2, for the PSCz (upper histogram) and 1.2 Jy (lower, shaded histogram) surveys. We define a selection function. 
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</>(r), as the fraction of galaxies that can be seen out to a redshift distance r = (expressed in h^^ Mpc), given 
the survey flux limit. To determine 4'{r) we use the parametric maximum likelihood estimate proposed by Yahil et 
al. (1991), in which the following analytic form for the selection function is assumed: 

0(r) = Ar~'"^ (^TT") ^ 

The value of the normalization constant, A, docs not affect the modeling of peculiar velocities. In this paper we 
arbitrarily set 0(r) = 1 if r < Ts with Tg = 6 h^^ Mpc. This choice is equivalent to imposing a lower cutoff in the 
60 luminosity which, in turn, avoids giving too much weight to faint, nearby IRAS galaxies that may not trace 
the galaxy distribution reliably in the nearby volume (Rowan Robinson et al. 1990). The relevant parameters have 
been determined via likelihood analysis using only the objects within 100 hT^ Mpc . The results are displayed in 
Table 1 for the PSCz and 1.2Jy surveys, both in redshift (z-) and real (R-) space (i.e. after correction for redshift 
space distortions as discussed in §3.1). The number density of galaxies as a function of redshift distance, predicted by 
eqn. (1), is shown in the lower part of Fig. 2, superimposed on the observed N{r) histograms. The theoretical curves 
provide a good description of the data. The ratio of the number of observed to the number of expected galaxies gives 
the radial overdensity field, shown in the upper part of Figure 2 for both the PSC^: (thick line) and the IRAS 1.2Jy 
catalogues. A horizontal line is drawn at the mean density. The radial density fields traced by the two surveys are 
consistent with one another, but the larger shot noise in the 1.2 Jy sample makes the amplitude of density fiuctuations 
larger than in the PSC2: survey at large radii. 

Our estimator for the selection function is independent of clustering but requires prior knowledge of the evolu- 
tionary rate. Springel and White (1998) have recently developed a new technique to estimate the rate of evolution 
of the selection function. For the 1.2Jy catalogue they find rather strong evolution and this is even more dramatic 
in the PSCz catalogue (Springel 1996). Ignoring this effect could introduce spurious streaming motions in the model 
velocity field. However, we have verified that our selection function for PSCz is very similar to the one derived by 
Springel (1996) within the scales relevant to our analysis. The difference between the two selection functions increases 
with distance but it is only 5% at 200 Mpc. Strong evolutionary effects become important only beyond these 
scales and we shall therefore ignore them in our modeling. 

For the Abell/ACO composite sample we adopt the selection function derived by BP96: 

1 if r < Tci 

0(0 = I 0.5(1 + ^1 exp(-r/roi)) if r,i < r < r,2 (2) 

0.5(^2 exp(— r/ro2) + A\ exp(— r/roi)) if r > Tci 

The estimated values of the parameters are listed in Table 2, and the expected number of clusters as a function of 
redshift distance is shown in the lower plot of Fig. 2 as a thick line. 

The moan separation, I, of the objects in a population of tracers limits the intrinsic resolution with which we 
the population probes the underlying cosmic fields. Thus, Rs = is a natural smoothing length which keeps the shot 
noise at a constant level. In a magnitude limited survey, I increases with distance r according to: 

Kr) = [0(r){n)]-^/^ (3) 

where the average number density of objects (n) has been estimated as 

1 

= (4) 

i=l 

and the sum extends over all Ng objects contained within the volume sampled, V . The value of the average number 
density within 100 hT^ Mpc , computed in redshift and real space, is shown in Table 1 both for the PSCz and IRAS 1.2 
Jy catalogues. The more physically meaningful values of the mean galaxy separation at 100 Mpc , l{r) expressed 
in Mpc~^, are also given in the table. Because of our normalization (0(r) = 1 within 6 Mpc), our estimate of 
(n) for the 1.2 Jy galaxies is slightly different from that of Fisher et al. (1995a). Fig. 3 shows the mean interobject 
separation as a function of r for the three samples considered. Because of the shallower depth of the 1.2Jy sample. 
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its intergalaxy separation increases much more steeply with distance than in the PSCz; sample. Thus, at a fixed 
resolution, the PSCz survey probes cosmic structures out to larger distances than the 1.2 Jy survey without an 
increase in shot noise. The dot-dashed line in Fig. 3 shows the Abell/ACO mean intercluster separation. Locally, this 
is much larger than the corresponding separation in the galaxy surveys, but on scales above 170 Mpc , clusters 
start to become more effective than IRAS 1.2Jy galaxies in tracing the cosmic density fields. 



3 NON-PARAMETRIC MODELS OF THE COSMIC VELOCITY FIELD 

The main aim of this paper is to obtain a self-consistent, non-parametric model of the velocity field in the local 
universe (i.e. for r < 150 Mpc). We do this by removing the redshift space distortions that affect the observed 
radial galaxy distribution in the survey using two different procedures. The first is the iterative technique devised 
by Strauss & Davis (1988), applied by Yahil et al. (1991), and further refined by Sigad et al. (1998). The results 
presented here were obtained by applying the original Yahil et al. (1991) technique. The second procedure, used here 
primarily as a check for possible systematic effects, is the non-iterative technique developed by Nusser and Davis 
(1994). 



3.1 The methods 

The two procedures that we have implemented rely on three important assumptions. Firstly, cosmic structures are 
assumed to have grown by gravitational amplification of small amplitude fiuctuations present in the density field at 
early times (gravitational instability, c.f. Peebles 1980). Secondly, fluctuations on the scales of interest are assumed 
to be small enough that linear (or mildly non-linear) theory is applicable. Thirdly, luminous objects (galaxies or 
clusters) are assumed to trace the underlying mass density field according to: 

Hr^)=b,5{v), (5) 

where 5 is the mass overdensity at position r, 5i is the fluctuation in the number of luminous objects at this same 
location, and hi is the biasing parameter. Eqn. (5) is often referred to as the "linear biasing model." The reconstruction 
methods return a model of the cosmic velocity field that depends on the parameter (3 — 0.%^ /hi. The value of (3 must 
be determined a posteriori by comparison with other observational data. 

The two reconstruction techniques that we use are extensively discussed and tested by Branchini et al. (1998) 
where the reader may find further details. Here, we simply outline the main principles of each method. 
Method 1: The iterative technique 

In linear theory the overdensity is related to the peculiar velocity field by 

V-v = ~l35i. (6) 
The model velocity fleld is then obtained by iteratively solving the system of equations: 

and 

r = cz — Y ■ {-v — -vlg)- (8) 

Eqn (^ follows from integrating eqn (|^; v is the peculiar velocity of a tracer object at position given by the radial 
unit vector f ; and vlg is the velocity of the Local Group. 

During each iteration, the selection function and the mean number density of the population of displaced objects 
are recomputed. Eqns (^) and (^) are only valid in the linear regime and so the force fleld needs to be smoothed 
to eliminate nonlinear effects. We employ a "top hat" fllter with a smoothing radius Rs{r) = max[5, {{n)(f>{r))^^^^] 
/i~^Mpc. This choice eliminates most of the nonlinear contributions and keeps the shot noise at a constant level. 
To improve numerical convergence, we adiabatically increase the value of (3 at each often iterations, from 0.1 to 1.0 
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(Strauss et al. 1992b). Note that the amphtude of the velocity vectors scales, to first order, linearly with /3. Unless 
otherwise stated, all the plots of the model velocity fields that we present assume /3 = 1. 

Around high density regions, such as clusters, there may be triple-valued regions in which the same redshift is 
observed at three different positions along a given line-of-sight. We correct for this by collapsing the galaxies within 
clusters and implementing the "robust procedure" ("method 2") of Yahil et al. (1991) to determine the locations of 
galaxies within triple-valued regions. For the collapsing procedure, we have used two different datasets: the six nearest 
clusters listed in Table 2 of Yahil et al. (1991) and a larger catalogue of 61 objects obtained by merging the 59 SMAC 
clusters of Hudson et al. (1998) with the Yahil et al. (1991) dataset. The results show that the collapsing procedure 
only affects the infall pattern around rich clusters. However, with the level of smoothing used here, the effect is 
negligible and no systematic differences are found when the model velocity fields are compared on a point -by-point 
basis. 

Method 2: Spherical harmonics expansion 

This procedure is based on the method proposed by Nusser and Davis (1994). In linear theory, the velocity field 
in redshift space is irrotational and thus may be derived from a velocity potential: 



Expanding the potential and the galaxy density field, 5", in spherical harmonics, and using the Zel'dovich approxi- 
mation, it can be seen that the two fields obey a Poisson-like equation: 



where 5f^ and are the spherical harmonic coefficients of the galaxy density and velocity potential fields respec- 
tively. We have considered the coefficients up to i < Imax = 15 and |m| < Imax. We employ a Gaussian filter with 
smoothing radius Rs{r) — max\^.2,{(ji)<j){r))~^^^] Mpc , compute the coefficients in redshift shells for the 
different catalogues, solve eqn. (^) for ^im, and then compute the full three-dimensional velocity field in redshift 
space using eqn. (^. Galaxy positions and peculiar velocities at the real space positions, r, are obtained by assuming 
a one-to-one mapping between the real and redshift space positions, r and s, along the line-of-sight. We minimize 
the problem of triple-valued regions by adopting the same cluster-collapsing scheme used above. 

Redshift space distortions are only one of several effects that hamper the recovery of cosmic density and velocity 
fields from observational data. Incomplete sky coverage is another, potentially serious, problem and we deal with 
it using a filling method similar to that introduced by Yahil et al. (1991). We use their original technique to fill in 
the region at galactic latitude |6| < 8°, but each real PSCz galaxy at \b\ < 8° is included to replace a synthetic 
object that resides in the same longitude-distance bin. Masked regions at larger galactic latitudes are filled in with a 
random distribution of synthetic galaxies having the observed mean number density. Tests performed by Branchini et 
al. (1998) using the mock catalogues described in the next section show that the only bias introduced by this filling 
in technique is a spurious bulk motion of < 60 kms~^(see §6.1). As a further check we have also used the new Fourier 
interpolation technique developed by Saunders et al. (1998b), which interpolates structure across the mask under the 
assumption that the density field of galaxies is a lognormal random field. In doing this, we have used a 'tapered', or 
radially-dependent, mask to take into account the known incompleteness in the PSCz at large distances (> 15000 
kms~^) and low latitudes {Ab > 1™), due to difficulties in making secure identifications and obtaining redshifts for 
obscured distant galaxies. The area of sky across which we interpolate is defined by = 2'" at V < 15000 km s~^, 
increasing to Ab = 1™ at 1/ = 30000 kms~^. Outside this volume, the PSCz is essentially complete. The model 
velocities predicted using this mask turned out to be almost identical to the ones obtained using the simpler filling 
in procedure described above. 

Another source of random and systematic error is the radial selection function which, when coupled to redshift 
distortions, generates the so called 'rocket effect' discussed by Kaiser and Lahav (1989). This is a spurious acceleration 
measured from a magnitude-limited sample of galaxies by an observer who has a peculiar velocity unrelated to the 
true gravitational acceleration. In Method 1 this effect is quantified and corrected for using the mock catalogues 
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discussed below. In Method 2 the rocket effect is explicitly accounted for by the second term in the right-hand side 
of eqn. (0). 

Finally, another potential source of bias is the fact that IRAS galaxies are preferentially of late morphological 
type and are less clustered than early types (Strauss et al. 1992a, Hermit et al. 1996). However, any biases arising 
from this are likely to be small since Baker et al. (1998) have shown that the ratio of the amplitudes of IRAS and 
optical galaxy density fields, smoothed on scales > 5/i~^ Mpc, is almost constant, implying a ratio of linear biasing 
parameters of ho/hi ~ 1.4. 

3.2 Error estimates using mock catalogues 

We have used a suite of large cosmological N-body simulations to generate mock galaxy catalogues that mimic the 
main properties of the PSCz and 1.2Jy redshift surveys. We use these to quantify random and systematic errors in 
our reconstructed velocity fields. For the fields inferred from cluster catalogues, we adopt the error estimates derived 
by BP96 and by Branchini et al. (1997) on the basis of a hybrid Monte Carlo/mock catalogue technique. 

The simulations we use are those described by Cole et al. (1998). We consider two different cold dark matter 
cosmologies: a flat model with S7o = 0.3 and cosmological constant term, Ac^/3Hq = 0.7, and a critical density 
universe {flo = 1.0) with power spectrum shape parameter, F = 0.25. In both cases, the fluctuation amplitude 
was normalized by the observed abundance of galaxy clusters. This requires setting as = O.SSr^^'^'® (Eke, Cole & 
Frenk 1996), where ag is the linear rms mass density fluctuation in top hat spheres of radius 8^~^ Mpc. Galaxies 
were identified with random particles from the simulations (so that 6 = 1 by construction). Prior to generating the 
mock catalogues, we smoothed the velocities in the simulations using a top hat filter of radius 1.5 ^~^Mpc. This 
brings cri2(l), the pairwise velocity dispersion of objects with projected separation < 1 Mpc , down to ^ 250 
kms~^(when only particles within 30 Mpc from the observer are considered). This value is in accordance with 
the recent analysis of the Optical Redshift Survey of Santiago et al. (1995, 1996) by Strauss, Ostriker and Cen (1998) 
for galaxies outside high density regions, and with the estimate of cri2(l) for late type galaxies in the Perseus-Pisces 
redshift catalogue by Guzzo et al. (1997). Both these measurements refer to 'field' galaxies, and are therefore well 
matched to the late type IRAS galaxies of the PSCz catalogue. 

Each mock PSCz catalogue is contained in a sphere of radius 170 Mpc and, for the purposes of error analyses, 
10 nearly independent mock catalogues were obtained from each cosmological model. Several constraints were applied 
in order to obtain mass distributions and velocity fields as similar as possible to those observed: 

1) Local Group observers. Hypothetical observers were selected from the set of particles with velocity, v_lg ~ 625 ±25 
kms~^, lying in regions in which the shear within 5 Mpc is smaller than 200 kms^^and the fractional overdensity 
within the same scale ranges between -0.2 and 1.0 (Brown and Peebles 1987). These constraints mimic the Local 
Group environment. 

2) Coherent galaxy dipole. A galactic coordinate system, {l,b), in the (periodic) computational volume was defined 
such that the velocity of the observer pointed towards {l,b) — (276, 30), the direction of the dipole anisotropy observed 
in the cosmic microwave background (CMB) radiation. This dipole is known to be approximately aligned with the 
dipole in the distribution of IRAS galaxies (Strauss et al. 1992a, Schmoldt et al. 1998). 

3) Radial selection. We generated flux-limited "galaxy" samples using a Monte Carlo rejection procedure to select 
particles in the simulations, assigning them fluxes according to eqn. (1). In the vicinity of an observer, the particle 
number density in the simulation is less than the number density of galaxies in the PSCz and 1.2Jy catalogues. We 
therefore generated catalogues that are semi-volume limited at a radius of 10.9 Mpc for PSC2: and 7.8 Mpc for 
1.2Jy. 

4) Masked areas. To mimic the incomplete sky coverage, we excluded all objects within the unobserved regions in the 
two IRAS catalogues. 

The distribution on the sky of the galaxies in one of our PSCz mock catalogues is shown, as an Aitoff projection, 
in the lower panel of Fig. 1, where it may be compared to the real survey displayed in the upper panel. An illustration 
of the utility of these mock catalogues is provided in Fig. 4. The upper left-hand panel displays a two dimensional slice. 
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corresponding to the mock Supergalactic plane, through the density and velocity fields within a region of 120 Mpc . 
Both fields were tabulated on a 64^ grid and subsequently smoothed using a Gaussian filter with a constant effective 
radius of 6 Ti"'^ Mpc . The upper right-hand panel shows density and velocity fields reconstructed using Method 1. All 
the main features, as well as most of the small scales structure, are correctly reproduced. The agreement between the 
fields can be better appreciated from the two lower plots. The one on the left shows the absolute value of the errors 
in the reconstructed density field. The main discrepancies occur within the zone-of-avoidance, depicted as a shaded 
region, and at large distances where the shot noise is large. A correlation between the error and the signal is present 
only in the outer regions. The lower right-hand plot displays the differences between the two velocity vector fields. 
Again, the largest discrepancies occur at large distances and no strong coherence is detected. Very similar results are 
obtained when the reconstruction is performed using Method 2. 

Unless otherwise stated, the same Gaussian filter of 6 /i^^ Mpc is used throughout the rest of this paper. As 
discussed in §2, a smoothing radius of 6 Mpc corresponds to the average galaxy separation at ~ 80 Mpc, i.e 
to the maximum distance out to which comparisons with observed peculiar velocities are still possible. The drawback of 
using a constant smoothing is that the shot noise increases with distance, artificially enhancing the density contrast 
in the outer regions, as may be seen in Fig. 4. With our adopted smoothing, the average shot noise in 5{r) at a 
distance of SO/i"^ Mpc is ^ 96%, rapidly increasing to ^ 180 % at 120 Mpc . A complete description of the mock 
catalogues, their use for error estimation and a detailed assessment of the reliability of our galaxy reconstruction 
methods is given in Branchini et al. (1998). The corresponding analysis for the Abell/ACO cluster catalogues may 
be found in Branchini et al. (1997). 

4 A COSMOGRAPHIC TOUR 

In this Section wo present a qualitative description of the model density and velocity fields of the local universe 
derived from the PSCz survey. We analyse the data using Method 1 above. Visualizing three dimensional structures 
and the corresponding vector fields, is not easy. Fortunately from this point of view, the most interesting structures in 
the nearby universe, on scales larger than the Local Supercluster, are roughly distributed along a planar structure, de 
Vaucoulcurs' (1953) Supergalactic plane (at Supergalactic coordinates SGZ= 0). In this analysis we will mainly follow 
the custom of displaying the density and velocity features in a two dimensional slice through this plane. However, 
distributions along other parallel planes will also be displayed. 

With the depth and sampling frequency of the PSCz survey, a reliable map of the density field can be constructed 
out to a distance of 150 Mpc . Fig. 5 shows a slice through the Supergalactic plane of am adaptively smoothed 
ovcrdcnsity field. Within 30 Mpc a constant Gaussian filter of radius 3 Mpc was used, but beyond that the 
radius of the filter increases linearly with distance up to 11.25/i~^ Mpc at 150 Mpc . This variable smoothing filter 
was chosen to maintain a roughly constant sampling error over the large range of scales considered. It represents a 
compromise between resolution in the inner regions and acceptable shot noise errors near the edge. The drawback is 
that, unlike the case of constant smoothing, the density contrast decreases with distance, as is clear from the map. 
(The S — level is indicated by the yellow line.) The most striking feature of this map is the large-scale coherence 
of the mass distribution. Interconnected overdensities, separated by very large voids, extend over distances of order 
100 Mpc . The most prominent structure, which plays a major role in the dynamics of the local flow fleld, is the 
overdense ridge extending from the Perseus-Pisces supercluster, close to the centre of the map, all the way to the 
Shapley concentration near the top left corner. 

The PSCz survey is large enough that the velocity field corresponding to the mass distribution can be recon- 
structed quite accurately, with a relative error always smaller than 50% in the region depicted in Fig. 5. This precision, 
however, is higher than that that of measured galaxy peculiar in a volume this size and so we cannot compare our 
model predictions with velocity data over the entire region of Fig. 5. A reliable comparison is only possible over a 
smaller volume, typically of radius ~80/i~^ Mpc . For this reason we now describe in some detail our reconstructions 
within this distance. 
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4.1 The PSCz model density and velocity fields within 80 h ^ Mpc 

Fig. 6 shows the PSC2: model density and velocity fields smoothed with a 6h^^ Mpc Gaussian in a slice through the 
Supergalactic plane. With this smoothing, the Local Supercluster, centered on the Virgo cluster at (SGX, SGY) = 
(—2.5, 11.5), does not appear as an isolated structure but is connected instead to the prominent Hydra-Centaurus 
supercluster at (SGX, SGY) = (—35,20). Together with the Pavo-Indus-Telescopium supercluster [(SGX, SGY) = 
(—40, —15)], the latter makes up the well known Great Attractor. The Coma cluster and its neighbour, A1367, appear 
as a peak at (SGX, SGY) = (0, 75), slightly elongated in the SGX direction. The Perseus Pisces supercluster, clearly 
visible at (SGX, SGY) = (45, —20), is the second largest peak on the map, well separated from its northern extension 
[(SGX, SGY) = (45, 20)] which is sometimes called the Camelopardalis supercluster. Finally, the Cetus Wall may be 
seen as an elongated structure around (SGX, SGY) = (15, -50). The Sculptor void [(SGX, SGY) = (-20, -45)] is the 
largest underdensity in the map, but is almost matched in size by a void in the background of the Camelopardalis 
supercluster. Throe more undcrdcnso regions that exert an important influence on the local dynamics arc the voids in 
the foreground of Coma [(SGX, SGY) = ( — 10, 50)], in the background of the Perseus-Pisces complex [(SGX, SGY) = 
(50, -50)], and behind the Great Attractor [(SGX, SGY) = (-60, 10)]. 

The competing dynamical roles of the various overdense and underdense structures seen in Fig 6. determine 
the local velocity field in a complex manner that cannot be simply described as a bulk flow or by a multi- spherical 
infall model. The local velocity field implied by the PSCz density field is illustrated by the vectors plotted in Fig. 6. 
Its dominant features are the large infall patterns towards the Great Attractor, Perseus-Pisces, Cetus Wall and 
Coma. A striking property is the large-scale coherence of the velocity field, apparent as a long ridge between Cetus 
and Perseus-Pisces and as a large-scale flow along the Camelopardalis, Virgo, Great Attractor baseline and beyond 
(sec §5). A prediction of the PSCz data is the lack of prominent back- infall onto the Great Attractor region. The 
flow around it appears to bo determined by the compressional push of two voids (at (SGX, SGY) = ( — 10,50) and 
(—60, 10)) and the pull of the Shapely Concentration on much larger scales (see §5). All these features are also present 
when Method 2 is used for the reconstruction (§4.2), or when the 1.2Jy survey is used as the input catalogue (§4.3). 
As pointed out by Davis, Nusser & WiUick (1996), this model velocity field does not match the Mark III peculiar 
velocity field (WiUick et al. 1997a) which exhibits a largo outflow away from the Centaurus supercluster and an inflow 
onto the Hydra complex. There is also no evidence for a motion of the Perseus Pisces supercluster towards us. 

Fig. 7 extends our qualitative analysis to two planes of constant SGZ, 40h~^ Mpc above and below the Super- 
galactic plane. The slice at SGZ= -|-40 Mpc shows a region which is largely underdense (Fig 7b), dominated by a 
void which is connected to the Local Void identified by Tully (1987). The prominent peak at (SGX, SGY) = (40, —30) 
seems to be an extension of the Perseus-Pisces supercluster and the one at (SGX, SGY) = (—50,20) appears to 
be connected to the Pavo-Indus-Telescopium complex. Overall the density fleld traced by the PSCz survey looks 
remarkably similar to the one obtained from the IRAS 1.2Jy catalogue (see §4.3 and the analogous discussion by 
Strauss and WiUick 1994). The velocity field (Fig 7a) is still characterized by a coherent, large-scale, flow towards 
the same direction [(SGX, SGY) = (—50, 50)] as the stream seen in the Supergalactic plane. At negative SGZ= —40 

Mpc the extensions of the Pavo-Indus-Telescopium (—50,15) and Perseus-Pisces (20,-20) superclusters are 
visible. The dynamical effect of these two peaks is evident in the associated large infall patterns (Fig. 7c). 



4.2 Comparison of the two reconstruction methods 

A detailed quantitative error analysis of mock PSCz catalogues performed by Branchini et al. (1998) also shows that 
the predicted velocity field is not affected by any significant systematic bias. To check this result further we have 
compared the reconstructions produced by the two different methods discussed in §3.1. The two upper plots in Fig. 8 
show the predicted velocity fields smoothed with a 6h~^ Mpc Gaussian in a slice through the Supergalactic plane. 
The two methods succeed in reproducing the main features that we have already highlighted: remarkable large-scale 
coherence in the velocity field, clear infall patterns onto Coma, Perseus-Pisces and the Great Attractor but no back- 
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infall onto the latter. The only noticeable difference is that Method 2 seems to blur slightly the sharp features (like 
the Cetus ridge) produced by Method 1 in the velocity field. 

Quantitative evidence for the similarity of the two model fields is provided in the lower part of Fig 8, which shows a 
scatterplot of the SGY Cartesian components of the reconstructed velocity fields, smoothed with a 6/i~^ Mpc Gaussian, 
measured within 80/i~^ Mpc. For both model velocity fields we have used (3 — 1. In the plots we show only 1000 
points taken at random from the 73490 available. The peculiar velocities in the two models should obey the linear 
equation, 

M2 = BMi + A, (11) 

where Mi and M2 denote any of the Cartesian components of the peculiar velocity predicted by Methods 1 and 2, 
respectively, and B is expected to be 1 if no systematic errors are present. The offset A allows for possible differences 
in the predicted bulk flows in the two models. We obtain the values of A and B by minimizing 

2 _ (Ma.. -A- BKh,if 

where omiA and (TM2,i are the errors in the velocities at a generic gridpoint i in the two methods, and Nt is the total 
number of points used for the comparison. We assume that ajv/^.i — o-M2,i = o", so that eqn.( |l^ ) becomes 



2 _ I (M2,. -A- BNh,,f _ xl 

^ _2 / ^ /-I I T>2\ _2 • \ "^1 



cr2 ^ (1 + 52 



Only Ni of the A^t points used in the comparison are independent. It can be shown that the quantity Xe// = {Ni/Nt)x^ 
is approximately distributed as x^ with Nd.o.f — Ni — 2 degrees of freedom (e.g. Hudson et al. 1995). We can therefore 
approximately evaluate a by setting Xeff= Nd.o.f so that 



^ iV.(iV.-2)- 

We take Ni to be the number of independent volumes within the volume sampled. With our smoothing radius, 
Ni ~ 6300. We then find a regression slope, B= 0.98 ± 0.06. The zero point, A= —47 ± 5, indicates a systematic 
difference in the SGY component of the bulk velocity in the two models. This offset is only evident when the two 
velocity fields are compared in the CMB frame and almost disappears in the Local Group (LG) frame. As explained 
in §3, the peculiar velocities are reconstructed in the LG frame. Transformation to this frame is performed by adding 
the reconstructed LG velocity at each galaxy position. This can lead to systematic differences if the LG velocities 
predicted by the two methods do not agree. Indeed, the zero point offset measured in the velocity-velocity scatterplot 
for the three Cartesian components, (-123,-47,132)/3 kms~^, is similar to the difference between the predicted LG 
velocities in the two models, (-150,-50, 98)/3 kms~^. No other systematic errors are detected as the regression slopes 
are close to one also for the two remaining components (0.97 ± 0.06 for SGX and 1.06 ± 0.05 for SGZ.) For the 
SGX and SGZ components, the dispersions around the fit are only slightly larger than for SGY {a = 61, 54, 70)/3 
kms~^respectively.) This represents the intrinsic error in the reconstruction and is nearly a factor of 2 smaller than the 
average total error derived from the error analysis on the mock PSCz catalogues (which also accounts for uncertainties 
in the filling in procedure, nonlinearity, finite volume effects, etc) by Branchini et al. (1998). 



4.3 Comparison of the PSCz and 1.2Jy model velocity fields 

Since they were drawn from the same parent catalogue, we expect the PSCz and 1.2Jy catalogues to give consistent 
model velocity fields, at least in the nearby volume where the sampling by 1.2Jy galaxies is not too sparse. In 
analogy with the PSCz fields displayed in Fig. 6, Fig. 9 shows the 1.2Jy density and velocity fields, smoothed with 
a 6/i~^ Mpc Gaussian, in a slice through the Supergalactic plane. Most of the characteristic structures previously 
identified in Fig. 6 are also visible in Fig. 9 except that underdense regions in the 1.2Jy map appear somewhat more 
extended than in the PSCz map. The overall pattern in the PSCz velocity field is reproduced in the 1.2Jy map 
although the large-scale coherent flow along the Camelopardalis, Great Attractor baseline is less evident in the latter. 
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On the other hand, the infall pattern around the Great Attractor in the 1.2Jy survey spreads out over a larger region. 
These discrepancies might well be due to the larger shot noise in the 1.2Jy catalogue. 

Our 1.2Jy model velocity field is consistent with that derived independently by Webster, Lahav and Fisher (1997) 
using the method developed by Fisher et al. (1995b). Comparison with their Fig. 7d (which assumes a CDM model 
with r = 0.2 as a prior in the Wiener filtering procedure) reveals only one noticeable difference between the two 
maps. This is in the region of the Great Attractor, where the model of Webster et al. predicts a weak back-infall. The 
discrepancy, however, is small and may simply refiect our use of a constant 6h~^ Mpc Gaussian smoothing, which is 
somewhat larger than the smoothing applied by Webster et al. 

The lower part of Fig 10 displays the difference between the PSC^; and the IRAS 1.2Jy overdensity fields, displayed 
in Figs. 6 and 9 respectively, A = Spscz — Si.2jy.. Positive (continuous lines) and negative (dashed lines) contours are 
shown, with the thick line delineating the A = contour. The main discrepancies between the two fields occur within 
the zone-of-avoidance, in the Great Attractor region and beyond. Outside the Galactic plane, the radially increasing 
shot noise is responsible for the large discrepancies seen in the outer regions. Apart from these local features, we 
notice that the mean value of A, (A), is positive in the region shown in the figure. Indeed, this is true within the entire 
spherical volume of radius 70 Mpc in which we find {Spscz) = 7 - 10~^ and (5i.2jy) = —6-10^^ /i^Mpc"^ (for both 
catalogues the average density has been measured within 100 /i^^Mpc). If the errors were due to sparse sampling 
alone, then this discrepancy would be significant at about the 2a level. However, the errors in the estimates of the 
average density are larger than this. In particular, as pointed out by Davis, Strauss and Yahil (1991), the intimate 
coupling between the estimation of the selection function and the average density itself generates uncertainties in the 
reconstructed average density which are large enough to account for the discrepancy we have detected. The mismatch 
in the average densities on a scale of 70 Mpc generates the infall pattern visible in the upper panel of Fig 10, 
while smaller scale features originate from local mismatches between the two density fields. 

Fig. 11 (lower panel) shows a more quantitative comparison between the PSCz and 1.2Jy overdensity fields, both 
smoothed with a 6h~^ Mpc Gaussian and tabulated onto 64^ grids. The comparison is made in the CMB frame and 
assumes j3 = 1. Only those gridpoints within 80 Mpc of the LG position are considered. As in the previous section, 
we fit a straight line to the data and estimate the parameters by minimizing . The slope of this line is 1.08 ± 0.8, 
indicating that there are no systematic differences between the two density fields. The negative zero point offset, 
—0.5 ± 0.03, reflects the discrepancy between the average densities of the two catalogues discussed in the preceding 
paragraph (A5 = —0.05). The scatter as = 0.20 is similar to the average shot noise at the gridpoint positions both for 
PSC^; {asN = 0.18) and for the IRAS 1.2 Jy {asN = 0.24) catalogues. Similar considerations apply to the velocity- 
velocity comparison in the upper panel in which, as in Fig 8, we only show the SGY Cartesian component of the 
velocity. The comparison, performed in the CMB frame, assumes /3 = 1. There is a non-negligible zero point offset of 
(56 ± 8)/? kms~^which is comparable to the difference between the SGY components of the predicted LG velocities 
for two catalogues (80/3 kms"'^). The dispersion around the fit, cr„ = 104/3 kms~^, is again similar to the average 
shot noise in the peculiar velocities, computed as in Yahil et al. (1991) (88/3 kms'^for PSCz and 140/3 kms'^for 
IRAS 1.2Jy). Similar results are found from the scatterplots of the two remaining Cartesian components. Their slopes 
are consistent with unity (1.06 ± 0.6 and 0.91 ± 0.08 for SGX and SGZ respectively) and the zero point offsets are 
comparable to the discrepancies arising from the transformation from LG to CMB frames (153/3 and 118/3 kms~^). 

Finally, Tadros et al. (1998) have recently suggested that the PSC^: catalogue should not be used down to 
0.6 Jy and that a more conservative cut in flux, 0.745 Jy, should be adopted instead. To check whether possible 
incompleteness at the low flux limit could affect the modeling of the PSC«: velocity field we have generated several 
different models for density and velocity fields using a suite of PSCz subcatalogues taken at different flux limits. With 
a 0.75 Jy cut the density and velocity flelds are still remarkably similar to those with the cut at 0.6Jy. Increasing 
the flux limit maJies the PSCz density and velocity fields more and more like the IRAS 1.2 Jy fields with which 
they almost coincide when the cut is taken at 1.2 Jy. Overall, our prediction for the PSC2: velocity and density fields 
are robust and are not affected by incompleteness in the catalogue at low flux levels, at least within a region of 200 
/i-^Mpc. 
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5 MODELLING VERY LARGE SCALE MOTIONS 

The sampling density of the PSC2; survey is high enough to allow investigation of the density field out to ^ 120 
Mpceven with a 6/i~^ Mpc Gaussian smoothing. This is illustrated in Fig. 12 which displays the usual section 
of the density and velocity model along the Supergalactic plane. The shot noise in this map is the same as in 
the maps discussed in §4.1. The larger volume reveals the full extent of the coherent streaming involving galaxies 
from the Camelopardalis Supercluster all the way to the Shapley Concentration (that begins to become visible at 
(SGX, SGY) = ( — 100, +60)), passing through the Local Supercluster and the Groat Attractor. 

The dynamical sources of the coherent motion seen in Fig. 12 are the same ones that we identified earlier as 
being responsible for the flow pattern behind the Great Attractor. The gravitational pull is mostly due to the Shapley 
Concentration, but the coherence is aided by the combined push of the two large voids (at (SGX, SGY) = (—80, 10) 
and (SGX, SGY) = ( — 50, 70)) that straddle the ridge and were only partially visible in Fig. 6. This map also reveals 
a largo extension to other undordensc regions, the Sculptor void and the two connected voids behind the Cetus- 
Perseus-Pisces-Camelopardalis complex. 

In Fig. 13 we show slices along four planes parallel to the Supergalactic plane, at SGZ= ±40 and SGZ= ±80. 
Fig. 14 shows the corresponding velocity fields. The same structures that we had identified in Fig. 7 can now be traced 
over larger scales. The shot noise is still small and so we can fully appreciate the large scale coherence of features in 
both the density and velocity fields. For example. Figs. 13a and 13b show that the void that we had identified in the 
upper half of Fig. 7b actually extends over a much larger volume. Similarly, the stream towards (SGX,SGY) = (-50,50) 
that we had identified in the SGZ=0 and SGZ=+40 Mpc planes persists at SGZ=+80 Mpc. No very large 
scale features are present in the plane at SGZ=-40 Mpc (Fig.l3c and 14c) which, as we had previously noted, 
is dominated by a few isolated infall patterns onto the southern extensions of the Pavo-Indus-Telcscopium and 
Perseus-Pisces superclusters. New large coherent structures, however, begin to appear again in the plane at SGZ=-80 
/i-^Mpc. 

5.1 Comparison with the Abell/ACO model fields 

Various surveys of IRAS and optical galaxies have been used to construct model velocity fields which proved to be 
remarkably consistent with one another (Yahil et al. 1991, Kaiser et al. 1991, Hudson 1994, Baker et al. 1998). We 
carry out a similar exercise here, using a completely different set of mass tracers, the Abell/ACO clusters. With these 
we can model the density and velocity fields out to very large scales, albeit with much larger sampling errors. BP96 
have already employed the same Abell/ACO cluster subsample described in §2 to model the density and velocity 
fields out to 250 Mpc . The reconstruction technique they used is a simplified version of our Method 1 (cf. §3.1) in 
which the selection function (which is nearly constant on the scales of interest) is not iteratively updated, no special 
treatment is given to triple-valued regions and the procedure for filling in the zone-of-avoidance (|6| < 20°) is the one 
used by Yahil et al. (1991). In order to compare the PSC^; and cluster model fields we smooth both with a Gaussian 
filter of width 20 Mpc . Such a large smoothing is required because of the large intercluster separation. 

Fig. 15 shows the smoothed density fields along the Supergalactic plane, within a distance of 150/i~^Mpc, 
derived from the PSCz (upper panel) and cluster (lower panel) catalogues. Overdensity contours are plotted in steps 
of AS = 0.2 for the galaxy field, and of AS = 0.88 for the cluster field. This is equivalent to rescaling the cluster field 
by a factor be = 4.4 which, as we shall see in §5.2 below, is our estimate for the relative linear bias parameter between 
IRAS galaxies and galaxy clusters. As before, the heavy line traces the S = contour. The dashed lines show the 
approximate location of the zone-of-avoidance in the PSC« (|&| < 8°) and cluster (|6| < 20°) samples. 

Despite the be = rescaling, the density peaks appear more prominent in the cluster sample than in the PSC« 
map. This may reflect more severe undersampling of the cores of rich clusters by IRAS galaxies than was indicated by 
Baker et al. (1998), or it may be due to shot noise in the cluster sample which, in spite of the heavy smoothing, is still 
substantial. It should be noted, however, that the Great Attractor region (SGX, SGY) = (—50, 0) and the Perseus- 
Pisces supercluster (SGX, SGY) = (50,-20), where the effect is strongest, both lie within the zone-of-avoidance of 
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the cluster sample. It is therefore possible that the height of these peaks has been artificially amplified by the coupling 
between shot noise and the filling-in procedure. For this reason we shall exclude the region at \b\ < 20° from the 
quantitative analysis in §5.2 below. Above |6| = 20° the Shapley Concentration (SGX, SGY) = (—120, 70) is the only 
peak with a larger amplitude in the cluster map than in the PSCz map. Amplitudes aside, the positions of the peaks 
in the two density maps are very similar. The low density regions are also approximately coincident, although they 
are less extended in the cluster map in which voids also appear somewhat shallower. The latter effect simply refiects 
the fact that the assumption of local biasing, Sc = bc6g, breaks down at low Sg because of the constraint 6c > —1. 

Similar considerations apply outside the Supergalactic plane. Fig 16 extends the comparison between clusters 
(left-hand side) and PSCz galaxies (right-hand side) to four slices parallel to the Supergalactic plane, at SGZ=±45 
fr^ Mpcand SGZ=±100 h'^ Mpc . Although some discrepancies in the position of peaks and voids exist, the large- 
scale features are remarkably similar across all the slices. In particular, the large overdense region at positive SGZ 
is evident in both the clusters and PSCz density fields. Unlike in the Supergalactic plane, there appears to be no 
appreciable systematic difference between the amplitudes of density peaks of clusters and IRAS galaxies. 

Fig. 17 shows the velocity fields inferred from the two samples, again with a 20h~^ Mpc Gaussian smoothing, in a 
slice along the Supergalactic plane. The amplitude of the velocity vectors is on an arbitrary scale but, for the clusters, 
the velocities have been rescaled by the factor be = 4.4. The two large overdensities lying within the zone-of-avoidance 
in the lower panel of Fig. 15 largely determine the infall patterns at those locations [(SGX, SGY) = (—50, 0)] and 
[(SGX, SGY) = (50,-20)]. Beyond |6| = 20°, however, the galaxy and cluster velocity fields are very similar. Both 
exhibit a large coherent flow along the Camelopardalis, Great Attractor, Shapley Concentration baseline. The infall 
onto the Shapley region is more prominent in the cluster map but, in general, the velocity field patterns are very 
similar at positive SGY, with an outfiow at positive SGX, and a convergent fiow towards the Shapley Concentration 
at negative SGX. Below the zone-of-avoidance, at negative SGY, both maps show outflow from the Sculptor void, 
but this is less prominent in the clusters map. 

Note that a similar good agreement would have been obtained by considering the density-velocity maps by 
Scaramella (1995) which also use the Abell/ACO cluster distribution. 

5.2 The relative linear bias between galaxy clusters and IRAS galaxies. 

If biasing is a local process then, in the regime where mass fluctuations are small, we expect the amplitude of the 
bias to be independent of scale (see e.g. Cole et al. 1998). This expectation is consistent with the results of POTENT 
analyses performed by Dekel et al. (1993) and Sigad et al. (1998) and we expect it also to be valid in the present 
analysis in which the PSCa and cluster density and velocity fields have been smoothed with a 20h~^ Mpc Gaussian 
filter. Thus, on the scales of interest, we expect the two model density and velocity fields to be linearly related: 



where C and P stand for cluster and PSC« and represent either S or any one of the Cartesian components of the 
velocity field. The constant be is the bias parameter of clusters relative to IRAS galaxies and Ac allows for a relative 
offset in the mean density or bulk velocities of the two fields. For quantitative analyses it is most convenient to use 

the SGY Cartesian component of the velocity field since this is the least affected by uncertainties in filling in the 
zone-of-avoidance. To estimate be we adopt the strategy of Hudson et al. (1995) and Branchini et al. (1997), already 
used in §4.2, of regressing the two model fields by minimising the quantity 



where the subscript i refers to any of the Nt gridpoints within the comparison volume. The quantities ac,i and ap^i 
represent the errors in the cluster and PSC2: fields, respectively. 

The errors in the cluster field, crc,i, have been estimated by Branchini et al. (1997). They are the sum in quadrature 
of the intrinsic errors in the reconstruction procedure, as estimated by BP96 using Montecarlo techniques, and the 



C = Pbc + A. 



(15) 




(16) 
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shot noise uncertainties which are dominant for the very sparse cluster fields. The latter were evaluated using the 
mock catalogue generated by Kolatt et al. (f996) which is designed to reproduce the distribution of structures in 
our local universe. A typical error in the cluster overdensity field is {(Jcs) = 0.36, and in the SGY-component of the 
velocity field it is (ac.vy) = 250/3 kms~^. 

The errors in the PSC2 fields have been estimated using the mock catalogues described in §3.2. The basic 
procedure consists of comparing the density and velocity fields obtained by applying Method 1 to the mock galaxy 
catalogues with the true fields in the parent N-body simulation. Fig. f 8 (which is analogous to Fig. 4) compares the 
true and reconstructed density and velocity fields of a mock PSCz catalogue (different from the one shown in Fig. 4) 
in a slice along the Supergalactic plane (upper panels). The maps of the density and velocity errors are shown in the 
lower panels. Both fields are smoothed with a Gaussian filter of radius 20 Mpc . As in Fig. 4, the errors in the 
density and velocity reconstruction occur preferentially within the zone-of-avoidance. However, the correlation with 
the signal, which was present in the outer regions of Fig. 4, has been almost erased by the large smoothing. We have 
characterized the reconstruction errors by noticing that the residuals correlate with distance, galactic latitude and, for 
the velocity field, with the signal itself. From the analysis of the mock catalogues, we have derived two approximate 
expressions for the errors (valid for > 20°): 

{ap,s) = O.f 1 - 6.7 fO""* ■ \b\ + 3.3 fO""* • r (17) 
and 

(crp,,„) = 70 + 2.7 fO"' ■ r + 0.f5 • \vy\, (f8) 

where r is the distance measured in h"^ Mpc , 6 is the galactic latitude in degrees (not the linear biasing factor), and 
\vy\ is the amplitude of the SGY Cartesian component in kms^^. The actual assumed error at a gridpoint has been 
generated from a Gaussian distribution centred on (crp) with dispersions of 0.04 for the overdensity and 40 kms~^for 
the Vy field, respectively. 

The regression of the PSCz and cluster fields uses all the gridpoints within f20 Mpc and outside the cluster 
zone-of-avoidance (i.e. at |ti| > 20°). As a result of the large smoothing applied, not all the gridpoints in the comparison 
volume are independent. The number of independent points, Ni, can be computed as in Dekel et al. (f993), 

Nt Nt 

= ^M-rl/2Rl), (f9) 

j=i i=i 

where Tij is the separation between gridpoints i and j and Rs is the smoothing radius of the Gaussian filter. This 
expression has been derived for a smooth density field and its application to a velocity field requires some caveats. 
For a given comparison volume, the large coherence of the velocity field causes the number of independent points to 
be smaller than A'^;. On the other hand, our peculiar velocities are predicted using a sample covering a larger volume 
than that used in the comparison and this increases the number of independent points. Because of these competing 
effects, we simply approximate the number of independent points by Ni, even in the v — v comparison. 

As in §4.2 we define the Xe// = {Ni/Nt)x'^ statistic which corresponds, in practice, to multiplying the errors ap 
and ac by Nt /Ni in equation (^6|). We use this statistic to estimate the errors in be and Ac- 

Fig. f9 (upper panel) shows a 5-S scatterplot of the model cluster and PSCz overdensity fields measured at 
~ 1000 randomly selected gridpoints. All 6426 original points are used in the regression analysis. The solid line 
shows the best fit obtained by minimizing Xe// ^^nd the parameters of the fit are listed in Table 3. The resulting 
bias parameter is 6f — 4.4 ± 0.6. This is consistent with results from an independent likelihood analysis in which the 
phenomenological power spectra of IRAS galaxies and Abell clusters were compared (e.g. Peacock and Dodds 1994). 
The parameter ~ Xeff l^dof = 1.14, which may be taken as an indication that the errors have not been grossly over- 
or underestimated. The systematic difference in the amplitude of the density peaks in the cluster and IRAS 5 fields, 
noticed in Fig 15, manifests itself as a deviation from the best fitting line at large Spscz- Restricting the regression to 
values of Spscz < 0.35 returns 6f = 4.5±0.6, almost identical to the previous value, and 5* = xij j / Ndof ~ 1.15. This 
suggests that the exact weighting of PSCz galaxies in high density regions has only a minor effect on the regression 
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analysis, mainly because the overdensity mismatch in large density peaks is restricted to a very few gridpoints. 
As a further check, we have repeated the 5-S regression using a PSCz velocity model derived without applying 
our standard procedures for collapsing clusters and handling triple valued regions (which should exacerbate any 
discrepancies associated with high peaks.) The results, listed in the second row of Table 3, show that the effect on be 
is indeed very small, leading to fef — 4.5 ± 0.6. 

The lower panel of Fig. 19 shows the scatterplot of the SGY Cartesian components of the two velocity fields in 
which, as in the upper panel, to avoid overcrowding, we only show ~ 1000 randomly selected gridpoints. As shown 
in Table 3 the slope of the best fitting line is fee" = 4.7 ± 0.6 with no significant zero point offset, indicating that 
the SGY components of the cluster and PSCz bulk flows are consistent with one another over the scales of interest. 
However, as indicated in Table 3, the resulting Xe// large, S"" — 1.55. This could be due to an underestimate 
of the errors since the error analysis by Branchini et al. (1997) is based on clusters from the Kolatt et al. (1996) 
mock catalogues and these do not accurately match the Abell/ACO cluster distribution and velocities. We can obtain 
5"" ~ 1.0 if we allow for a reasonable error underestimate in the cluster field of ^ 30%, in which case we obtain 
be" = 4.0 ± 0.6. Note that consistent values of and 6c" are obtained within 1-a whether or not the cluster errors 
have been underestimated at this level. 

It is worth emphasizing that the agreement of the 5-5 and Vy-Vy comparisons is not trivial. The 5-5 comparison 
is local; it is hardly affected by problems related to filling in masked regions but is potentially prone to the cluster 
core weighting problem. The Vy-Vy comparison, on the other hand, involves the distribution of objects within the 
entire sample and thus is much more strongly affected by the unknown mass distribution within the zone- of- avoidance 
and beyond the sample's edge. We might therefore expect the two comparisons to be subject to different biases. The 
agreement in the estimates of be from the two analyses suggests that systematic biases have been properly taken 
into account and that the linear biasing assumption is a good approximation, at least on scales larger than our 20 

Mpc smoothing. 



6 THE BULK VELOCITY VECTOR 

In this section we consider the peculiar velocity, Vi,{R), of spheres of radius R centred on the Local Group. This is 
a low order statistic that, in principle, can be estimated observationally. The expectation value of the bulk velocity, 
\Vb\, averaged over scale R is: 

i\MR)f) = ^ J P{k)W{kRfdk, (20) 

where P{k) is the power spectrum of density fluctuations, W{kR) = 3 ain(kR)-kRcos(kR) ^j^^ Fourier transform 
of the spherical top hat window function in real space, chosen to facilitate comparison between the theoretical 
definition (eqn. and its observational analogue (eqn. ^). It is worth stressing that eqn.( assumes equal 
volume weighting while the use of different weighting schemes may lead to quite different results (Kaiser 1988, 
Strauss et al. 1995, Giovanelli et al. 1998). If the initial fluctuation field obeys Gaussian statistics, then evolution 
through gravitational instability preserves a Gaussian distribution for the amplitude of each Cartesian component 
of Vb{R), so that |Wi,(i?)| has a Maxwellian distribution. This property makes it difficult to constrain P{k) from the 
measured Vb{R). Nevertheless, comparison of the measured Vb{R) with the velocities predicted from the PSCz gravity 
field allows, in principle, an estimate of the (3 parameter. In practice, however, the bulk velocity is extremely sensitive 
to systematic errors both in the observational data and in the models. We attempt to take this carefully into account 
in the following analysis. 



6.1 The model bulk flow 

The basis for our treatment of random and systematic errors in the model bulk velocity vector are, again, the mock 
PSCz catalogues described in §3.2. Finite volume effects are not a concern here since the volume of the simulation 
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box is comparable to that of the PSCa sample. However, neglecting modes on scales larger than the size of the sample 
can bias the comparison between the model bulk flow and the measured one. In what follows we will treat these two 
effects separately. We first account for the intrinsic errors in the bulk fiow model and we then address the problem of 
minimizing the effect of missing large-scale power. 

From each mock catalogue we generate a model velocity field using Methods 1 and 2 of §3.2. The field is then 
smoothed onto a 64^ cubic grid of side 192 ft ~ ^ Mpc using a 12ft~^ Mpc Gaussian filter. The smoothing filter used 
here is different from the one used before because our aim is to perform a homogeneous comparison with the Mark 
III and SFI bulk flows. We measure the cumulative bulk velocity vector in the CMB frame by averaging over the 
peculiar velocities measured at gridpoints: 

= ^i^,'^)<n^^'^'\ (21) 

where Wij^fc is the predicted velocity vector in the CMB frame at gridpoint {i,j,k). The sum X](ijfe)<fl extends 
over all the gridpoints contained within a sphere of radius R. The same exercise is repeated using the original N- 
body velocity field. This gives an unbiased estimate of the true bulk velocity and the comparison between true and 
reconstructed velocities is an estimate of the error in the PSCz model bulk flow. 

In the upper panel of Fig. 20 we show the average difference between reconstructed and true cumulative bulk 
velocity vectors. Each Cartesian component is displayed with a different symbol. The errorbars, representing the dis- 
persion around the mean, are shown only for the SGY component. The dispersions for the remaining two components 
are larger by nearly a factor of two. From the plot we see that the errors in the SGY and SGZ components of the 
reconstructed bulk flow are mainly random errors, but the SGX component is affected by a systematic positive bias. 
The amplitude of this systematic error decreases with distance (from ~ 61 kms~^at 10 ft~^Mpcto ~ 16 kms~^at 
100 Mpc ). Very similar results are obtained if the reconstruction is carried out using Method 2. These systematic 
errors arise mainly from the fiUing-in procedure for the zoric-of-avoidance which, as wo discussed in §3.2, affects both 
the velocity fields and, via uncertainties in modeling the LG velocity, the transformation from LG to CMB frame. To 
corroborate this hypothesis we have performed 10 reconstructions using unmasked mock catalogues. With full sky 
coverage, the systematic bias in the SGX component disappears and the dispersions in the SGX and SGZ components 
decrease by a factor of two. The systematic error in the SGX component is ultimately induced by the requirement of 
having only LG-like observers in the PSCz mock catalogues. In particular, visual inspection reveals that the galaxy 
dipole constraint in §3.2 almost invariably implies the existence of a Great Attractor-like density peak at a distance 
of 30-60 Mpc , lying at negative SGX and partially overlapping the zone-of-avoidance. The fiUing-in procedure 
seems to underestimate the amplitude of this partially hidden density peak and this in turn causes a systematic offset 
in the SGX component of the reconstructed peculiar velocities, with minor effects on the two other components, as 
shown in Fig. 20. A similax bias is to be expected in the true model PSCz bulk fiow because the Great Attractor 
really exists in our universe. To obtain bias-free estimates of the PSCz bulk flows, we must correct for such systematic 
errors. Clearly, the best way of circumventing the problem is by restricting attention to the SGY component of the 
bulk velocity. However, sometimes one is interested in the amplitude of the bulk velocity and, in this case, a direct 
numerical correction of the systematic error is appropriate. 

Since the amplitude of the estimated errors depends on /3, the errors displayed in the upper panel of Fig. 20 (or 
the analogous quantity for the velocity amplitude) cannot be directly used to make quantitative corrections to the 
PSCz model bulk flow. However, we can take advantage of the fact that, to flrst approximation, the model velocities 
scale linearly with /3. Thus, the ratio between reconstructed and true bulk velocities should be independent of /3. 
This ratio, averaged over our mock catalogues, is shown by filled circles in the lower panel of Fig. 20, as a function 
of distance, with errorbars again representing the dispersion around the mean. Since it is independent of (3, this ratio 
may be used as a multiplicative factor to correct the predicted PSCz bulk flow amplitudes for systematic errors. As 
expected, the corrections become very small when considering the SGY component (open circles). We have performed 
a similar analysis to estimate errors in the PSCz bulk velocity obtained using Method 2 and in the 1.2Jy bulk flow 
model, obtained using Method 1. 
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We are now ready to derive model bulk flows and their uncertainties from the PSCz and 1.2Jy surveys. The 
upper left panel in Fig 21 shows the amplitude of the cumulative bulk velocity vector, predicted using various samples 
and methods. In all cases the bulk flow was computed by integrating the density field out to the same distance of 200 
hr^ Mpc , hence enforcing the same volume bias in all cases. The velocities are normalized to /3 = 1. Open and filled 
circles show PSCz results using Methods 1 and 2 respectively, and filled squares show results from the 1.2Jy survey, 
all corrected for systematic errors as discussed in the preceding paragraph. For clarity only the 1-cr error bars from 
Method 1 are plotted. The filled triangles show the bulk velocity computed from the Abell/ACO model velocity field, 
rescaled by he = 4.4 (see §5.2). There is remarkably good agreement between the cumulative bulk flows computed 
from the different sets of mass tracers, with typical deviations of less than ^ 20% from the mean in the amplitude 
from the different surveys and analysis methods. The directions of the cumulative bulk ffow vectors, measured within 
60 Mpc, are plotted in the upper right-hand panel of Fig. 21. The PSC2; model bulk ffow points, within la, to 
the direction of both the 1.2Jy and Abell/ACO model bulk flows. The direction of the CMB dipole, indicated by the 
asterisk, is plotted for reference at the centre of the figure. In the lower panel of Fig. 21 we show the misalignment 
angle of the various bulk fiows with respect to the PSCz -Method 1 model. In all cases the misalignment is small 
and almost independent of radius. The only exception is the increasing misalignment angle at large radii between the 
PSC2 bulk velocities determined using Method 1 and Method 2. Note the close alignment between the PSCz and 
clusters bulk fiows at all radii. 

So far we have modelled bulk velocities in the nearby universe by considering only the gravitational effect of the 
mass distribution within 200/i~^ Mpc . Bulk ffows, however, are sensitive to the mass distribution on scales larger than 
those probed by our samples. Neglecting large scale modes leads to a volume bias which can affect the comparison 
with observations. Various methods have been proposed to restore the missing large scale modes (e.g. Strauss et al. 
1995, Tormen and Bertschinger 1996, Cole 1997). Here we propose a simple statistical treatment based on linear 
theory, using eqn. ( po| ) to obtain a correction factor for the missing contribution to the amplitude of the bulk velocity 
from scales beyond the sample boundary. The correction factor is the ratio 

P{k)W{kRfdk 

F{P{k),k^,^) = ^^^^^ p(^k)WikRrdk- ^^^^ 

The numerator in this equation is the mean true bulk ffow, while the denominator is the mean bulk flow generated 
only by density fluctuations on scales smaller than Rmax = In our case, Rmax ~ 200/i~^Mpc. Note that 

the ratio depends only on the spectral shape and on kmin', the dependence on /3 cancels out. Tests using large N- 
body simulations (e.g. Jenkins et al. 1998) have shown that this approach is indeed effective. Since the amplitude 
of the bulk ffow follows a Maxwellian distribution, and its components are Gaussian random variables, we can 
numerically estimate the uncertainty in this correction. Assuming that these uncertainties and the intrinsic random 
errors computed above are uncorrelated, we obtain the total random uncertainties on the model bulk flow by adding 
them in quadrature. In computing F we assume a CDM model with spectral shape F = 0.25, as suggested by a wide 
variety of data on the large-scale galaxy distribution (e.g. Baugh 1996). 

We have applied the correction of eqn. (22) to the bulk velocity reconstructed from the PSCz survey using 
Method 1. The resulting cumulative bulk velocity (for /3 = 1) is shown in Fig. 21 as a dot-dashed line labelled MIL. 
We regard this as our best estimate of the /3-dependent bulk flow. 



6.2 Model vs. observed bulk flow 

In this Section, we compare our best estimate of the predicted bulk flow velocity from the PSCz survey with recent 
observational estimates. This comparison serves two purposes. Firstly, consistency between predicted and observed 
velocities lends support to the hypothesis that structure grew by gravitational instability and gives confldence in the 
integrity of the observational data. Secondly, the comparison allows an estimate of the parameter /3 = Q,q'^ /b. 

The determination of bulk flows from peculiar velocity surveys is prone to systematic errors. For example, 
zero-point errors in the calibration of the distance indicators coupled with limited sky coverage may mimic a bulk 
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flow. Similarly, the coupling of large intrinsic errors in distance measurements with inhomogeneities in the density 
distribution (the inhomogeneous Malmquist bias) also results in a spurious outflow. Thus, to perform an unbiased 
comparison with theoretical predictions requires full-sky, homogeneous surveys of peculiar velocities and an accurate 
model of the survey's window function. 

In Fig. 22 we plot observational determinations of bulk flows derived from four (almost) independent datasets. 
The lower panel shows the full cumulative bulk flow as a function of distance while the upper panel shows its SGY 
component only. The filled triangles (taken from Dekel et al. 1998) represent the cumulative bulk flows in the Mark III 
catalogue. The sparse and noisy Mark III velocities have been smoothed assuming that the velocity field is irrotational. 
This guarantees that a unique three-dimensional velocity field is derivable from the observed radial velocities, as in 
the POTENT method (Dekel, Bertschinger & Faber 1990). The resulting three-dimensional peculiar velocity field, 
smoothed with a 12h~^ Mpc Gaussian, and defined on a regular grid is directly comparable to our model prediction. 
The filled square shows the preliminary result by Eldar et al. (1998) who use a POTENT smoothing technique to 
derive the bulk velocity from the SFI catalogue (Giovanelli et al. 1997a, 1997b). The open triangle in the lower panel 
is the bulk velocity inferred from 44 supernovae Type la by Riess, Press & Kirshner (1995), as reported by Dekel 
(1997) (the effective radius of this last measurement is much smaller than the depth of the sample because the data 
were weighted by the inverse of the errors). Finally, the open square at large distance shows the bulk velocity derived 
by Lauer & Postman (1994; LP94) from a sample of brightest cluster galaxies. The directions of the observed bulk 
velocity vectors are given in the right-hand panel. For the Mark III, SFI and Sn la determinations, they have been 
estimated from data within 50 Mpc . The direction of the LP94 bulk flow is taken from Strauss (1997) and refers 
to a depth of ~ 90 /t"^ Mpc. 

Except for the LP94 result, there is excellent agreement between the various determinations of the bulk flow in 
Fig. 22. These may be compared with the bulk flows predicted by the PSCz survey, indicated by the thick dot-dashed 
lines which enclose the 1-cr allowed range. These are normalized to /3 = 0.75 and corrected for volume bias. The filled 
circle in the right-hand plot marks the direction of the predicted bulk velocity vector within 50 Mpc . 

Requiring that the predicted bulk flow should match the measured one gives an estimate of f3. In performing this 
comparison we shall ignore the discrepant LP94 result which Strauss et al. (1995) and Watkins and Feldman (1995), 
amongst others, have argued is inconsistent with currently acceptable cosmological models. The LP94 data point 
is also inconsistent with our predicted bulk flow derived only from the gravitational instability and linear biasing 
hypotheses. Indeed, as shown in Fig. 22, matching the amplitude of the LP94 bulk flow would require a value of 
/3 ~ 1.9, which is incompatible with all other current measurements (e.g. Giovanelli et al. 1998a, 1998b, Dekel et al. 
1998 and Eldar et al. 1998) and implies a very large misalignment angle of 70° between the LP94 vector and the 
direction of our model flow velocity. In what follows we limit our comparison of model and observed bulk flows by 
imposing a series of restrictions designed to minimize possible systematic errors. These restrictions are: 

- We consider only the cumulative bulk flow from Mark III as estimated by Dekel et al. (1998). Like the PSCa model 
velocity field, the Mark III velocities have been smoothed onto a regular grid and filtered on a similar scale. This 
ensures that the comparison is as homogeneous as possible. 

- We consider Mark III peculiar velocities calibrated according to the recent VELMOD analysis by Willick and 
Strauss (1998). 

- We exclude scales smaller than 30/i~^ Mpc since small differences in the smoothing procedures applied to model 
and observed velocities in the nearby region can bias the comparison. 

- We do not include the result of Eldar et al. (1998) which, although consistent with those of Dekel et al. (1998), is 
still preliminary. 

We perform the comparison between model and observed bulk flows using a likelihood technique similar to that 
used by Strauss et al. (1992a) and Schmoldt et al. (1998). The aim is to estimate the likelihood of a particular value of 
(3 given the bulk flow observed in the Mark III catalogue and the one derived from the PSCa survey using Method 1. 
Under the assumption that the bulk velocities axe Gaussian random fields it is easy to show that the joint probability 
distribution for the observed and model bulk flows is a multivariate Gaussian that depends only on its covariance 
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matrix. The latter is completely specified by the value of 13 and by the power spectrum of density fluctuations, P{k). 
In this work we restrict attention to the family of CDM models for which P{k), parametrised as in Davis et al. (1985), 
is completely specified by its shape parameter, F = Q,oh, and normalization, erg. Although the latter can be regarded 
as an independent parameter, adopting the normalization inferred by Eke, Cole, and Frenk (1996) from the observed 
abundance of galaxy clusters, allows one to relate as, to Slo- In the following we consider only two cases, as, = 0.52 
and 0.87, which correspond to the cases of a critical, fio = 1.0, and an open, S7o = 0.3, universe. (For a flat Universe 
with S7o = 0.3, the required value of as is similar). 

Given a cosmological model specified by erg, /3, and F, we compute the joint probability of obtaining the cumulative 
bulk flow measured by Dekel et al. (1998) at 30, 40, 50 and 60 Mpc (vM(j), i = 1,4), and the PSCz; -Ml cumulative 
bulk flow displayed in Fig. 22, at all radii from 10 to 80 Mpc {yp{j),j = 1, 8). We ignore the correction for volume 
bias. The resulting probability density is: 

/(vM(i),i = l,4,vp(i),i = 1,8) = (27r)-'<''+**)/'(defAf)-'/"exp(-iv, •v™(M-i),™)dv„dvp, (23) 
where the covariance matrix Mim is: 



Mi™ = i < V, ■ v„ / dkP{k)Wi{k)W„,{k). (24) 

3 



The window functions are: 

Ti/ (i\ nSin(feri) - fcricos(fcri) 

WM(k) = 3 (25) 

for the Mark III cumulative bulk flow measured within n and 

~ sin(fcrj) -fcrjcos(fcrj) sin(fc_R,) - kRs cosjkRs) . , 

^^^^^-^ Ik^s 3 IkRj^ (26) 

for the PSCz cumulative bulk flow measured within rj. In eqn. (26), Rs = 200 Mpc represents the outer radius 
of the PSCi sample, i.e. the maximum distance out to which the PSCz density field is integrated to predict peculiar 
velocities. The second term in the right-hand side of eqn. (26) accounts for the volume bias affecting the Ml-PSCz 
bulk fiow. This correction is robust in the sense that the alternative strategy of limiting the integration in eqn. ( p4| ) 
to fc > fcmm ~ 2n/Rs gives very similar results. Finally, we take into account random errors in the measured and 
modeled bulk fiows by adding them in quadrature to the diagonal terms of the covariance matrix. 

From the joint probability density distribution (eq. |23| ) we construct the relative likelihood, C, of different world 
models (/?, F, ag): 

£=-21n(/). (27) 

The two plots in Fig. 23 show likelihood contours of £{vM{i),i = l,4;vp(j), j = 1,8) in the {l3,leq = l/F) plane, 
plotted at confidence levels of 68, 90 and 99 %. The shape of these contours is very similar for the two normalizations 
although the range of allowed values of /3 is somewhat larger for the high normalization. In both cases the likelihood 
peaks at /9 = 0.6 — 0.7; low values of /3 are excluded at high confidence level primarily by the requirement to match 
the amplitude of the observed bulk fiow. No stringent constraint can be imposed on the shape parameter F apart 
from an upper limit that follows from the fact that the power cannot be confined to very small scales. The cluster 
normalization also defines a (/3, F) {b, h) transformation. Therefore, we can explore the effect of introducing an 
observational constraint on the Hubble parameter, (0.4 < h < 0.8), which defines a vertical region in the {l3, F) plane. 
Introducing this extra constraint has a major effect in the case of low normalization since it forces the power to 
originate from unacceptably small scales. By contrast, in the case of high normalization, the extra constraint on h is 
perfectly consistent with the likelihood contours. 

Following Schmoldt et al. (1998) we can set a more stringent constraint on f} by considering not only the bulk 
flow but the entire set of observational information available on peculiar velocities in the local universe. In particular, 
we use the CMB dipole (vlg), the amplitude of the local velocity shear (|Vi,|), and the dipole in the gravitational 
acceleration induced by the distribution of PSCz galaxies in TV independent radial shells {vd{k),k — 1,N) in order 
to compute the constrained probability function f{vM{i),i = l,4;vp(j),j — l,8;Vtj(fc),fc = 1, N;-vlg\'Vs). The 
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modeling of this probability distribution, its covariance matrix and the treatment of the errors may be found in 
Schmoldt et al. (1998) and Strauss et al. (1992a) to which we refer the reader for a more detailed discussion. In this 
analysis we use the PSC2: differential dipole, Vd, measured by Schmoldt et al. (1998) in N = 15 non-overlapping 
top-hat windows of width 10 Mpcand the same values of vlg = 625 ± 25 and |vs| < 200 kms~^used in §3.2. 

The effect of adding these new observational constraints is shown in Fig. 24. The overall shape of the likelihood 
contours does not change appreciably. However, the valley along the F direction is now deeper and so the value of /3 
is better constrained. This is clearly seen in Fig. 25 where we plot the distribution of the normalized likelihood of /3 
obtained by marginalising the joint likelihood distribution over F. The continuous lines give the marginal distributions 
after integrating over all values of F, while the dot-dashod linos give the result of limiting F to the values allowed by 
the extra constraints on h. For ag — 0.87, we find (3 = 0.61q !J? (l^<^)i irrespective of the constraint on h. Thus, in a 
low density universe, a power spectrum normalized according to the cluster abundance is consistent with the velocity 
data. By contrast, for the low normalisation, as = 0.52, the constraints on /3 change appreciably with the addition 
of the new constraints, leading to a larger value of /? = 0.91° 20 (l-c)- As may be seen in Fig. 24, however, in this 
latter case the relative probability within the strip allowed by the constraint on is smaller than in the case of high 
normalization, thus rendering a high Qq universe less likely according to the velocity data. Finally, it is interesting 
to note that the assumption of linear bias coupled with the cluster normalization of the power spectrum implies 
that (3 ~ Q.5/ai^'~'^ , irrespective of the value of fio- Prom the variance of the PSC« galaxy counts at 8 Mpc, 
ffPSCz ^ Q ij^ obtain /3 ~ 0.7, consistent with the estimates above at the 1-cr level. 



7 DISCUSSION AND CONCLUSIONS 

The recently completed PSCa redshift survey of IRAS galax;ies represents an almost ideal dataset for studying the 
mass distribution and the gravity field in the local universe. In this paper we have used the PSCz survey to develop 
a nonparametric model of the local cosmic velocity field which supersedes results derived from the shallower 1.2 Jy 
and the sparser QDOT surveys. Our reconstructions are based on the assumptions that cosmic structure has grown 
by gravitational instability and that fluctuations in the galaxy distribution are proportional to fluctuations in the 
mass distribution. Wc have paid particular attention to a careful estimation of random and systematic errors using 
a suite of mock PSCz catalogues constructed from largo cosmological N-body simulations. As a further check on the 
validity of our results, we have implemented two independent methods for reconstructing the PSC« model velocity 
flelds, both of which give consistent results. 

Because of the large size of the PSCa survey, the density and velocity fields cam be reliably reconstructed out 
to a depth of 150ft~^ Mpc . Wo have presented maps of the galaxy distribution which, with a smoothing of XS 
(!)h~^ Mpc , clearly show the relative sizes of the main structures that characterize our local universe. The two largest 
peaks in our neighbourhood are the Great Attractor, made up of the Hydra-Centaurus and Pavo-Indus-Telescopium 
superclusters, and the Perseus-Pisces supercluster located in roughly opposite directions along the Supergalactic 
plane. The Local, ComarA1367, and Camelopardalis superclusters as well as the Cetus Wall are clearly visible in our 
maps, as is the giant Shapley concentration which appears near the edge of the survey, behind the Great Attractor. 
The laxgest underdensity in our vicinity is the well-known Sculptor void, but this is almost matched in size by a 
void in the background of the Camelopardalis supercluster. Three more underdense regions that exert an important 
influence on the local dynamics are the voids in the foreground of Coma, in the background of the Perseus-Pisces 
complex, and behind the Great Attractor. 

The local velocity field implied by the PSCz density field is complex. The dominant features are the organized 
infall patterns towards the large mass concentrations in the Great Attractor, Perseus-Pisces and Coma regions. 
Superimposed upon these axe impressive coherent flows along a ridge between Cetus and Perseus-Pisces and along the 
Camelopardalis-Virgo-Great Attractor-Shapley direction. We see no prominent back-infall onto the Great Attractor. 
Instead, the flow in this region is a result of an interplay between the compressional push of two nearby voids and 
the pull of the Shapley concentration on much larger scales. 
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The PSCa reconstruction of the velocity field agrees well with results from the 1.2Jy survey within 80h~^ Mpc , 
the region in which the latter provides adequate sampling. The only noticeable difference are the predicted bulk 
velocity vectors which (in the CMB frame) differ in the two surveys by ~ 130/3 kms~'^in each Cartesian component. 
This discrepancy arises from uncertainties in the way in which the zone-of-avoidance is filled in and from the sampling 
noise of the density field at large distances. While the former affects the PSCz and 1.2Jy results equally, the latter 
is more severe in the 1.2Jy case because of the larger shot noise. Thus, the misalignment between the Local Group 
acceleration and the CMB dipole vectors is ~ 25° for the 1.2Jy survey, but only ~ 15° for the PSC« survey (Schmoldt 
et al. 1998). Comparison of the spherical harmonic multipoles of the two model velocity fields in radial shells confirms 
the consistency of the PSCz and 1.2Jy survey predictions in the region of overlap (Teodoro et al. 1998). 

As was the case with the 1.2Jy survey, the velocity field reconstructed from the PSCz survey is partially incon- 
sistent with the peculiar velocities in the Mark III catalogue. A visual inspection of the PSCz and 1.2 ,]y velocity 
maps (Figs. 6 and 9) shows some difi^crcnccs around the Hydra Ccntaurus complex, where the Mark III data exhibit 
a large shear, and in the Perseus Pisces region, which, in the Mark III catalogue, appears to participate in a large 
streaming motion (Davis, Nusser and Willick 1996). As pointed out by these authors, beyond 30 Mpc there is a 
large coherent residual dipole in the difference between the 1.2Jy gravity field and the Mark III peculiar velocity field. 
This cannot be ascribed to improper modelling of the gravity field by IRAS galaxies since the dipole depends only 
on the mass within the surveyed volume. Furthermore, Baker et al. (1998) have shown that the gravity field in the 
optical ORS survey (Santiago et al. 1995, 1996) is consistent with the 1.2Jy gravity field (and thus, by implication, 
with the PSCz field as well). This suggests that the undersampling of cluster cores by the predominantly spiral galaxy 
population in IRAS surveys does not have a large effect in the inferred gravity field. The most plausible explanation 
for the discrepancy between the Mark III velocities and the IRAS (and ORS) model predictions seems to bo some 
systematic bias that affects the Mark III catalogue beyond 30 h^^ Mpc (Baker et al. 1998, Willick and Strauss 1998). 

As we saw in §5.1, the density and velocity fields inferred from the PSCz survey agree well with those inferred 
from a sample of Abell/ACO clusters. This is perhaps surprising since clusters are selected in a very different way 
from galaxies, but it is reassuring and suggests that systematic errors are under control in both cases. Comparison 
of the two model density and velocity fields, smoothed on the same cubic grid, out to a distance of 140 /i~^Mpc, 
gives the relative linear bias, foe, between the rich cluster population and PSCz galaxies. A simple analysis of 
the density density and velocity velocity comparisons gives very similar results, foe = 4.4 ± 0.6 and foe = 4.7 ± 0.6, 
respectively. These values are consistent with the estimate, foe = 4.5, obtained by Peacock and Dodds (1994) from an 
independent likelihood analysis of phenomenological power spectra. A more detailed comparison of the relative bias 
field of clusters and galaxies will be presented by Plionis et al. (1998). 

Finally, averaging over the peculiar velocity field reconstructed from the PSCz survey, we have calculated the 
expected bulk velocity of concentric spheres around us. Comparison of these predicted bulk flows with those measured 

from the Mark III catalogue by Dekel et al. (1998) gives an estimate of f3 — VP'^ /h. Wo implemented a likelihood 
analysis to carry out this comparison, taking into account the observed CMB dipole, the observed local shear field and 
the velocity of the Local Group predicted by the PSCz gravity field. If a high normalization for the power spectrum 
is assumed {as = 0.87), then our best estimate is /3 = 0.6lQ;jg (l-f). This value of (3 is consistent with results 
from analyses of the PSCz dipole (Schmoldt et al. 1998, Rowan-Robinson et al. 1998), and with the most recent 
determinations of P from velocity-velocity comparisons (Davis, Nusser & Willick 1996, Willick et al. 1997b, da Costa 
et al. 1998, Willick and Strauss 1998 and, within 1-a, also with Sigad et al. 1998). Adopting a low normalisation 
for the power spectrum gives a higher value of /3 = 0.91q Jo (l~o")- While still consistent with most of the estimates 
above, such large values are only marginally consistent with the measured variance of the IRAS galaxy density field 
and the power spectrum normalization derived by Eke, Cole and Prenk (1996). Furthermore, with this normalization 
the bulk flows arc also difficult to reconcile with the observationally favored values of Bubble's constant, /i, and the 
spectral shape paxameter, F. 

A more accurate determination of /?, to 15% accuracy, is possible by performing a more detailed comparison 
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between the peculiar velocity field inferred from the PSCa survey and the measured peculiar velocities at independent 
locations. An analysis of this kind is currently in progress. 
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Table 1. Parameters for the selection function of IRAS Galaxies 



Sample 


a 


/3 






/(lOO) 


{ 


n) 




z-IRAS PSCz 


0.54 


1.83 


6.0 


87.00 


11.67 


5.76 


10- 


-2 


R-IRAS PSCz 


0.52 


1.92 


6.0 


90.75 


11.66 


4.61 


10- 


"2 


z-IRAS 1.2Jy 


0.49 


1.80 


6.0 


51.27 


16.21 


6.41 


10- 


-2 


R-IRAS 1.2Jy 


0.47 


1.87 


6.0 


52.50 


16.24 


5.51 


10- 


"2 



Table 2. Parameters for the selection function of Abell/ACO clusters 



Sample roi 


ro2 Tci rc2 Al 


A2 


(n) 


Abell/ACO 31.8 


44.0 180 235 125 


289 


4.61 10-5 



Table 3. The clustcr/IRAS galaxy relative bias parameter, be- The top and bottom rows give results with and without applying 
the cluster collapsing procedure. Column 1: , the number of gridpoints used for the regression analysis; column 2: N^^^ f , 
the number of independent volumes; column 3: from the 5-S regression and it 1-(T error; column 4: A^, the zero point offset in 

the 5-5 regression and its l-cr error; column 5: = X^jj l^dof from the 5-<5 regression; column 6: 6^ from the Vy-Vy regression 
and its l-cr error; column 7: A^y , the zero point offset in the Vy-Vy regression and its l-cr error; column 8; 5"" = X^ff l^dof 
from the Vy-Vy regression. 





Nd.o.f. bi A« 


5* 6^" 


A"v 


S'"y 


6427 
6527 


54 4.4 ±0.6 0.06 ±0.07 

55 4.5 ±0.6 0.12 ±0.07 


1.14 4.7 ±0.6 

1.15 4.6 ±0.6 


-16 ± 90 
8.2 ±91 


1.55 
1.57 
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FIGURE CAPTIONS 

Figure 1. AitofF projection of the galaxy distribution in galactic coordinates, in the PSCz survey (upper panel) 
and in a mock catalogue constructed from a cosmological cold dark matter N-body simulation (lower panel). The 
flUed-in regions show unobserved or obscured regions; the zone-of-avoidance is the quasi-horizontal strip surrounding 
the galactic plane. 

Figure 2. Lower Plot: the number of galaxies as a function of redshift-distance in the PSCa (upper histogram) and 
1.2 Jy (lower, shaded histogram) samples. The curves show the expected counts as a function of distance estimated 
from the selection functions. The heavy line at the bottom shows the predicted distance distribution of Abell/ACO 
clusters. The labels give the total number of objects in each sample. Upper plot: the ratio between the observed and 
the expected number of galaxies in the PSCz catalogue (thick line) and in the IRAS 1.2 Jy catalogue (thin line). 
Figure 3. Mean inter-object separation as a function of radial distance (in Mpc ) in the three samples considered 
in this paper: PSCz galaxies (continuous line), 1.2 Jy galaxies (dashed line) and Abell/ACO clusters (dot-dashed 
line). 

Figure 4. Two dimensional slices through the density and velocity fields in a mock PSC« catalogue, corresponding 
to the mock Supergalactic plane. The hypothetical observer is located in a region analogous to the Local Group. 
Data are shown in a sphere of radius 12Qh~^ Mpc centred on the observer. Both density and velocity fields have been 
smoothed using a Gaussian filter of radius 6 h^^ Mpc . The upper left-hand plot shows the true fields while the upper 
plot on the right shows the reconstructed fields using Method 1. Continuous lines represent isodensity contours with 
a spacing of 0.5 in 5. Solid lines encompass overdensities and dashed lines underdensities. The heavy line indicates 
the (5 = contour level. The amplitude of the velocity vectors is on an arbitrary scale. In the lower left-hand plot we 
show the error map for the density field. Contours for the absolute value of the discrepancy between overdensities, 
A, are drawn at steps of 0.25. The shaded region around SGY=0 is the mock zone-of-avoidance. The lower plot on 
the right shows the difference between the true and reconstructed velocity fields. 

Figure 5. Real space density field derived from the PSCz survey. A slice along the Supergalactic plane is shown. 
The field has been smoothed with a variable Gaussian filter. The smoothing length is set at a constant value of 
3/i~^ Mpc within 30/i~^Mpcand increases linearly with distance up to a value of 11.257i~^ Mpcat 150 /i~^Mpc, 

where the most distant structure are located. The yellow line shows the 5 = contour. 

Figure 6. Real space density and velocity fields derived from the PSCz survey. The fields, smoothed with with a 
67i~^ Mpc Gaussian are shown in a slice along the Supergalactic plane. The most distant structures are located at 
80/i~^ Mpc . The thick continuous line shows the 5 = contour. Positive (continuous lines) and negative (dashed 
lines) contours are plotted at steps of A5 = 0.5. The amplitude of the velocity vectors, obtained for /3 = 1, is on an 

arbitrary scale. This reconstruction has been performed using Method 1. 

Figure 7. The same fields as in Fig. 6, but in two slices parallel to the Supergalactic plane. The top panels refer to 
the slice at SGZ= -|-40 Mpc , above the Supergalactic plane, while the ones at the bottom refer to the slice at 
SGZ= —40 Mpc , below the galactic plane. Density and velocity fields are plotted separately following the same 
conventions adopted in Fig. 6. This reconstruction has been performed using Method 1 for a value of /3 = 1. 
Figure 8. The peculiar velocity field reconstructed using Method 1 (top left panel) and Method 2 (top right panel), 
smoothed with a 6 Mpc Gaussian in a slice along the Supergalactic plane. The amplitude of the velocity vectors is 
on the same arbitrary scale in the two panels. The lower panel is a scatter plot of the SGY components of the velocity 
vectors illustrated in the upper panels. Only 1000 points, randomly selected from the ones within 80 Mpc, have 
been plotted. The parameters of the linear fit (A,B) are indicated in the legend, along with the scatter in the model 
velocities [a). 

Figure 9. The same as Fig. 6 but the density and velocity fields have been derived from the 1.2Jy survey. This 
reconstruction has been performed using Method 1. 

Figure 10. Map of the difference between the density and velocity fields obtained from the PSCz and the 1.2 Jy 
catalogues in a slice along the Supergalactic plane. The lower panel shows the difference between the two overdensity 
fields, A = dpscz — Si.2jy.: positive (continuous lines) and negative (dashed lines) isodensity contours are drawn at 
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steps of AA = 0.25; a thick line connects the points in which A = 0. The difference between the two velocity fields, 
both obtained for /? = 1, is displayed in the upper panel. The amplitude of the velocity vectors is on an arbitrary 
scale. 

Figure 11. Point-by-point comparison of the 1.2Jy and PSCz density and velocity fields within 80 Mpc . Only 
the values at 1000 randomly chosen gridpoints are plotted. The upper panel shows a velocity-velocity scatter plot 
for the SGY components of the peculiar velocity (assuming (5 = 1). The lower panel shows the 5-5 comparison. The 
parameters of the linear regression fit are given in both panels. 

Figure 12. A section, along the Supergalactic plane, through the real space density and velocity fields derived from 
the PSCz survey. The fields have been smoothed with a 6/i~^ Mpc Gaussian. This figure is similar to Fig. 6, except 
that it displays data in a larger spherical volume of radius 120/i~^ Mpc . The thick continuous line is the 5 = contour. 
Positive (continuous lines) and negative (dashed lines) contours are plotted at steps of A(5 = 0.5. The amplitude of 
the velocity vectors is on an arbitrary scale. This reconstruction has been performed using Method 1. 
Figure 13. Four slices, parallel to the Supergalactic plane, through the same density model illustrated in Fig. 12. 
The different panels refer to slices at SGZ=-|-40, -f-SO, -40 and -80 h^^ Mpc, respectively. The same conventions as 
in Fig. 6 are adopted. 

Figure 14. PSCz velocity fields corresponding to the density fields of Fig. 13. The amplitude of the velocity vectors 
is on an arbitrary scale. 

Figure 15. Real space density fields derived from the PSCz survey (upper panel) and Abell/ACO clusters (lower 
panel) within a distance of 140 Mpc . Both fields have been smoothed with a 20^~^ Mpc Gaussian. Here we show 
a slice along the Supergalactic plane. The 5 = level is indicated by the thick line. Other contours are plotted in steps 
of A(5 = 0.2 for the PSCz map and A5 = 0.88 for the clusters map. The dashed lines delineate the zone-of-avoidance 
in each sample. 

Figure 16. Cluster (left) and PSCz (right) density fields, as in Fig. 15, but in four slices along the Supergalactic 
plane, at SGZ=-f 100, +45,-45,-100 Mpc . The level of the density contours is also as in Figure 15. 
Figure 17. Peculiar velocity fields derived from the PSCz survey (upper panel) and Abell/ACO clusters (lower 
panel) corresponding to the density fields in Fig 15. Both fields have been smoothed with a 20/i~^ Mpc Gaussian. The 
amplitude of the velocity vectors is on an arbitrary scale, but the cluster field has been scaled according to a relative 
bias parameter, be = 4.4. The dashed lines delineate the zone-of-avoidance in each sample. 

Figure 18. Error maps similar to the one displayed in Figure 4, except that a different mock PSCz catalogue is used 
and a larger Gaussian smoothing of 20 Mpc is assumed. The density contours are drawn at steps of AA = 0.1 in 
the upper plots and of AA = 0.05 in the lower left-hand plot. 

Figure 19. Density and velocity scatterplots for the reconstructions based on the PSCz survey and a sample of 
Abell/ACO clusters. For clarity only ~ 1000 out of 1878 gridpoints within a sphere of radius 120/i~^ Mpc , and lying 
within |6| > 20°, are plotted. The velocity-velocity comparison in the upper panel refers to the SGY-component. 
The lower panel shows the 5-5 comparison. The parameters in the legend refer to the fits discussed in the text 
(see also Table |^. The errorbars give the mean 1-cr errors in the two model fields. 

Figure 20. Random and systematic errors in the model bulk flow derived from the analysis of mock PSCz catalogues. 
The upper panel shows the difference between the true cumulative bulk velocity vector and the vector reconstructed 
using Method 1. Different symbols are used for the three Cartesian components, as indicated in the figure. The lower 
panel shows the ratio between the amplitudes of the reconstructed and true bulk fiow vectors (filled circles) and their 
SGY Cartesian components (open circles). 

Figure 21. The amplitude (upper left-hand plot) and the direction (upper right-hand plot) of the cumulative bulk 
velocities predicted in different models. The circles give results from the PSCz survey using Method 1 (open circles) 
and Method 2 (filled circles); the filled squares give results from the 1.2Jy survey using Method 1; the triangles 
give results from the Abell/ACO cluster sample. For the models based on IRAS galaxies, the amplitudes have been 
normalized to /3 = 1, while for the cluster model the velocities have been normalized assuming a relative bias, be = 4.4. 
The errorbars give 1-a uncertainties obtained from the PSCz survey using Method 1. The dot-dashed line shows the 
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PSC« reconstruction using Method 1, statistically corrected to minimise the volume bias. The right-hand panel gives 
the direction of the predicted cumulative bulk flow at 60 Mpcin galactic coordinates {l,h). The asterisk at the 
centre marks the direction of the CMB dipole and the different samples are indicated by the same symbols used in the 
top left-hand plot. Errorbars are l-cr uncertainties derived from the analysis of the mock catalogues. The contours 
are set at constant misalignment angle from the apex of the CMB dipole, in steps of AO = 10°. The lower plot shows 
the misalignment angle between the various model bulk flows as a function of distance. The difference angle between 
models Ml and M2 (continuous line), Ml and 1.2 Jy (short dashed line) and Ml and Clusters (short dashed line) is 
displayed. 

Figure 22. The amplitude (left) and direction (right) of the cumulative bulk velocity vector. The bottom-left panel 
refers to the total velocity and the top-left panel to the SGY-component only. The filled triangles show the bulk 
velocity measured from the Mark III catalogue by Dekel et al. (1998) and the filled square the bulk velocity measured 
from the SFI catalogue by Eldar et al. (1998). The open triangle shows the bulk flow inferred from a survey of 44 
Type la supernovae by Riess, Press & Kirshner (1995). The open square gives the bulk velocity derived by Lauer 
& Postman (1994) from a survey of brightest cluster galaxies. The dot-dashed lines bracket the l-cr range of the 
bulk velocity predicted from the PSCz gravity field using Method 1, and corrected for the effects of long- wavelength 
modes. In the both panels the model predictions are normalized to /3 = 0.75. In the right-hand panel, the asterisk 
marks the direction of the CMB dipole in galactic coordinates and the flUed circle the direction of the bulk velocity 
at 50h~^ Mpc obtained from the PSC^; survey using the corrected Method 1. The other symbols correspond to those 
in the left-hand panel, at a distance of 50 hr^ Mpcin the case of the Mark III catalogue. 

Figure 23. Likelihood contours of C(-VM(i),i ~ l,4;vp(ji),j — 1,8) in the /3, leg — 1/T plane. The minimum (or 
the minima) is displayed by a small horizontal line. The left-hand panel assumes as = 0.87 and the right-hand panel 
(78 = 0.52. The vertical lines encompass the range for which 0.4 < h < 0.8, for f2 = 1 in the left panel and f2 = 0.3 in 
the right panel. 

Figure 24. As Fig. 23 but the contours displayed here refer to the conditional likelihood function £(vm(«),« = 
l,4;vp(i),j = l,8;Vd(A:),fc = 1, iV; v^gIv,). 

Figure 25. Marginal distributions of the likelihood of (3 for erg = 0.87 (left panel) and erg = 0.52 (right panel). 
Continuous lines give the result of integrating over all values of F, while dot dashed lines give the results of limiting 
the integration to the range encompassed by the vertical lines in Fig. 24. 
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